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Deep optical imaging of AGB circumstellar envelopes 
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ABSTRACT 

We report results of a program to image the extended circumstellar envelopes of asymptotic giant branch (AGB) stars 
in dust-scattered Galactic light. The goal is to characterize the shapes of the envelopes to probe the mass- loss geometry 
and the presence of hidden binary companions. The observations consist of deep optical imaging of 22 AGB stars with 
high mass loss rates: 16 with the ESO 3.6 m NTT telescope, and the remainder with other telescopes. The circumstellar 
envelopes are detected in 15 objects, with mass loss rates <; 2 x 10 -6 Mq yr -1 . The surface brightness of the envelopes 
shows a strong decrease with Galactic radius, which indicates a steep radial gradient in the interstellar radiation field. 
The envelopes range from circular to elliptical in shape, and we characterize them by the ellipticity (E = major/minor 
axis) of iso-intensity contours. We find that ~50% of the envelopes are close to circular with E 1.1, and others are 
more elliptical with ~20% with E ?Z 1.2. We interpret the shapes in terms of populations of single stars and binaries 
whose envelopes are flattened by a companion. The distribution of E is qualitatively consistent with expectations 
based on population synthesis models of binary AGB stars. We also find that ~50% of the sample exhibit small-scale, 
elongated features in the central regions. We interpret these as the escape of light from the central star through polar 
holes, which are also likely produced by companions. Our observations of envelope flattening and polar holes point 
to a hidden population of binary companions within the circumstellar envelopes of AGB stars. These companions are 
expected to play an important role in the transition to post-AGB stars and the formation of planetary nebulae. 

Key words. Stars: AGB and post-AGB — Stars: mass loss — Stars: circumstellar matter 



1. Introduction 

It is well established that stars lose a significant fraction of 
their mass during evolution on the asymptotic giant branch 
(AGB). The mass loss forms expanding envelopes of molec- 
ular gas and dust around the stars, and eventually recycles 
the stellar material back into the interstellar medium. The 
circumstellar envelopes are important for the study of many 
astrophysical processes, including circumstellar chemistry 
and dust formation. They are also important for stellar and 
galactic evolution, especially in connection with the forma- 
tion of planetary nebulae. For reviews on AGB envelopes 
and their properties see, e.g., Habing & Olofsson (2004) and 
Kerschbaum et al. (|201ip . 

Despite extensive study, there are many aspects of AGB 
envelopes that are not well understood. Here we focus on 
one of the most basic properties, the geometry of the en- 
velopes. On large size scales the morphology is known to be 
affected by interaction with the interstellar medium (e.g. 
Cox et al. 2012), but the shaping mechanisms close to the 
star are not well understood. It is not known if the mass 
loss of single stars is closely spherically symmetric, or if it 
is influenced by factors such as the stellar magnetic field. 
Binary companions are expected to influence the geometry, 
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Observatory, Chile (programs 078.D-0102, 082.D-0338 and 
0.84.D-0302) and on de-archived data obtained with the ESO 
Very Large Telescope and with the NASA/ESA Hubble Space 
Telescope 



but information on the presence of companions is problem- 
atic. Apart from symbiotic systems, only a few dozen AGB 
binaries are known (Jorissen I2008p . even though they are 
expected to be common. The reason so few are detected is 
because of the luminosity and variability of the AGB stars, 
and the opacity of their circumstellar envelopes. There must 
be a large population of hidden companions that can inter- 
act with the envelopes. 

The effects that stellar companions can have on AGB 
envelopes is suggested by hydrodynamical simulations 
(Theuns & Jorisse n [19931 Mastrodemos & Morris IT9991 
Gawryszczak 2002). There arc two main interactions. First, 
the reflex motion of the AGB star forms a spiral pattern 
in the envelope, and this has been detected in the case 
of AFGL 3068 (Mauron & Huggins[l!l hereafter MH06; 
Morris et al. I2006|) . Second, the large scale geometry of the 
envelope is expected to be flattened by a companion, and 
this has been discussed by Huggins et al. ([21)09) as a pos- 
sible probe of the presence of binary companions. Indirect 
shaping effects on the envelopes by jets launched from ac- 
cretions disks around a companion are also possible, as seen 
in the case of Mira (Josselin et al. ,2000). 

There is a lack of detailed observational information on 
the structure of circumstellar envelopes because the rela- 
tively cool material in the envelopes is difficult to image at 
high resolution. A review of techniques is given by Marengo 
(2009). Recent advances include millimeter interferometry 
of the molecular gas (e.g., Castro-Carrizo et al. I2010|) and 
imaging of the thermal dust emission, in the mid-infrared 
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Table 1. Properties of the AGB stars 



IRAS 


other name 


I 


b 


Chem. 


d 


y 


M 


/12 


f 12 / 125 


Refs. 












(pc) 


(kms -1 ) 


(10~ 5 M^yr -1 ) 


(Jy) 






00193-4033 


AFGL 5017 


326 


-75 





360 


17.3 


0.10 


312 


2.07 


1 


01037+1219 


IRC+10011 


128 


-50 





740 


21.6 


2.2 


1155 


1.19 


2 


04575+1251 


AFGL 5134 


188 


-18 





1200 


12.9 


1.0 


133 


1.09 


3, 4, 9 


05418-3224 




237 


-27 


C 


4300 


10.0 


1.1 


63 


2.47 


5 


06012+0726 


AFGL 865 


201 


-07 


C 


1470 


16.6 


1.8 


320 


1.42 


5 


06255-4928 




258 


-24 





1100 


15.0 


0.31 


66 


1.53 


4, 9, 10 


07454-7112 


AFGL 4078 


283 


-21 


c 


940 


13.5 


1.6 


613 


1.99 


5 


09116-2439 


AFGL 5254 


253 


+16 


c 


940 


12.8 


1.3 


737 


1.85 


5 


09452+1330 


IRC+10216 


221 


+45 


c 


120 


14.7 


1.75 


47500 


2.06 


6 


09521-7508 


AFGL 4098 


292 


-16 


c 


1300 


12.8 


1.1 


345 


1.86 


5 


10323-4611 




280 


+10 





580 


20.4 


1.0 


537 


1.46 


4, 9 


12384-4536 




301 


+19 





1100 


8.1 


0.20 


159 


2.25 


1 


14591-4438 




326 


+ 12 





600 


15.2 


0.25 


304 


1.47 


1 


16029-3041 


AFGL 1822 


345 


+16 





1900 


16.8 


1.10 


141 


0.53 


1 


17319-6234 


OH 329.8-15.8 


330 


-16 





1200 


17.5 


0.84 


252 


0.85 


1 


18240+2326 


AFGL 2155 


52 


+16 


c 


920 


15.1 


0.46 


731 


1.63 


5 


18467-4802 


OH348. 2-19.7 


348 


-20 





1200 


12.3 


0.78 


285 


0.83 


1 


19178-2620 


AFGL 2370 


12 


-18 





2200 


18.0 


1.3 


78 


0.81 


1 


20042-4241 


AFGL 5578 


358 


-31 





460 


14.4 


0.15 


221 


1.41 


7, 8 


20077-0625 


AFGL 2514 


36 


-20 





760 


16.0 


1.7 


1255 


1.18 


2 


23257+1038 


AFGL 3099 


92 


-47 


c 


1610 


10.5 


0.85 


190 


1.34 


5 


23320+4316 


AFGL 3116 


108 


-17 


c 


780 


14.7 


0.70 


959 


2.04 


5 



Refs: 

(1) d , V, and M from Loup et al. (1993 ) 

(2) d, V, and M from Olivier et al. (pOOTjl 

(3) V from OH observations, David et al. (1993) 

(4) M estimated through a correlation of M with the K— [12] color, where K is from the 2MASS catalog, and [12] = —2.5 
log(/i2/28.3), with /12 from the IRAS point-source catalog 

(5) d , V, and M from Groenewegen et al. I|2002p 

(6) d from Ramstedt et al. (PODS]) , M from Schoier & Olofsson (^D0T|| . Kx P from Loup et al. (TIMI)) 

(7) V from Loup et al. lfl993|l 

(8) d and M from Whitelock et al. (fl994)l 

(9) d estimated assuming L = 3500 Lq (e.g., Jackson et al. I2002p 

(10) No V cxp available, 15 kms -1 adopted value 



with ground-based observations (e.g., Lagadec et al. I201ip 
and in the far infrared with Herschel (Decin et al. 1201 ip . 
In this paper, we report progress using a different tech- 
nique: observing the envelopes in dust-scattered Galactic 
light. This technique makes use of deep optical imaging, so 
that high resolution is more readily attained than at longer 
wavelengths. 

Our initial observations using deep optical imaging pro- 
duced the first wide field image of the circumstellar enve- 
lope of the AGB archetype IRC+10216 at arcsecond reso- 
lution (Mauron & Huggins \MM 110001) • The core of the en- 
velope is illuminated by light from the central star, but the 
extended envelope is illuminated by the ambient Galactic 
radiation field, which reveals the detailed structure. We 
subsequently carried out pilot observations of other, more 
distant AGB stars, and found that the technique can pro- 
vide useful information on the shapes of the circumstel- 
lar envelopes (MH06). In this paper, we report more ex- 
tensive observations using the New Technology Telescope 
(NTT) at the European Southern Observatory at La Silla 
(Chile), with arcsecond resolution and exposures up to sev- 
eral hours. The goal is to extend the sample of AGB stars 
with resolved envelopes to determine their geometry. 



Section 2 describes the observations; Sect. 3 presents 
the images and information on the shapes, together with 
those of our earlier observations in the same format; and 
Sect. 4 discusses the results and their implications. 



2. Observations 

2.1. The sample 

The deep imaging technique is well suited for investigat- 
ing the shapes of the circumstellar envelopes of AGB stars 
with high mass loss rates. The thickest envelopes produce 
the maximum surface brightness and the largest envelope 
sizes for observation in dust scattered Galactic radiation. 
The high opacity of these envelopes also reduces or elimi- 
nates the stellar light contribution, which may be affected 
by structure close to the star and can mask the faint, ex- 
ternal illumination. 

The AGB targets for observations were selected from 
the catalog of Loup et al. (|1993j) . on the basis of strong 
IRAS fluxes and large infrared excesses (small /12//25 
ratios), as well as declinations accessible from La Silla. 
Preference was given to stars at relatively high galactic lat- 
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Table 2. Details of NTT observations Table 3. Observations from earlier work 



3 



TR 


Tnotr a 


XDcLIlLl C6 Eja L) US Lll (J 




Tl 






(min) 


("1 
\ ) 




nniQS— 4033 


II 


R (2W) V (150"l T (B) 


1.9 


(VI 


0457^4-1 ocjl 


II 


R fl 50^ V (QVl) T ("1 ") 


1.2 


(V) 


n^l 8— S994 


II 


R ^fiO'l J Ofta\ 

J-J 1 UU f 5 J. 1 Z/Uo 1 


2.2 


(R1 
y D ) 


0601 9-1-0796 


II 


T/ f49^ T (]()) 


1.6 


(VI 




II 


R (7u\ T (2u^\ 

i-J 1 1 OK ± I iUo 1 


2.2 


(R) 


07454—7112 


II 
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I 


R (301 /" (6l c 
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1.0 


(R) 


09521-7508 


I 


B (30) 


1.1 


(B) 


10323-4611 


II 


B (120), V (30) 


1.0 


(B) 


12384-4536 


II 


B (270), V (30) 


1.1 


(B) 


14591-4438 


II 


B (90) 


1.1 


(B) 


16029-3041 


II 


B (30), V (180) 


1.2 


(V) 


17319-6234 


II 


B (60), V (30) 


1.0 


(V) 


19178-2620 


II 


B (240), V (90), / (10s) 


1.7 


(B) 


20042-4241 


I 


B (60) 


1.4 


(B) 


23257+1038 


II 


V (510), 7 (6) 


1.6 


(V) 



a I = EMMI, II = EFOSC2 

6 FWHM of field stars in the images 

c /-band image obtained with EFOSC2 



itude, in order to minimize the effects of interstellar extinc- 
tion and crowded star fields. 

Details of all the AGB stars discussed in the paper are 
given in Table 1. The list includes new targets observed 
with the NTT (Section 2.2) and earlier observations which 
we include for completeness (Section 2.3). Table 1 gives the 
name, the Galactic co-ordinates, the chemistry (oxygen-rich 
or carbon-rich), distance d, expansion velocity T4 xp , mass- 
loss rate M, and relevant IRAS fluxes from the point source 
catalog. It can be seen that they represent AGB stars in the 
solar neighborhood (within a few kpc) with high mass loss 
rates, comparable in many cases with IRC+10216. 

2.2. NTT observations 

The observations with the NTT were carried out in four ob- 
serving runs. Two were made with the EMMI instrument 
in imaging mode on 12-14 October 2006 and 16-18 March 
2007, and two were made with the EFOSC2 instrument on 
20-23 February 2009 and 15-19 October 2009. The angular 
field is 5.'7 x 5'.7 with a pixel size of 0'.'3608 ± 0'.'0003 for 
EMMI, and 4.'1 x 4fl with a pixel size of 0'.'2406 ± 0'.'0002 
for EFOSC2. The first two observing runs were affected 
by poor weather, so most of the data reported are with 
EFOSC2, which also has greater sensitivity and better sam- 
pling of the PSF. 

Details of the observations including the instrument, the 
filters, and the exposure times are listed in Table 2. In many 
cases, we first secured an J-band exposure to ascertain the 
precise position of the target. The deep B or V images to 
detect the faint envelopes were made with a series of in- 
dividual exposures, typically of 30 min each, with a small 
shift in a or 8 between them to improve the quality of the 
final image. The data were reduced using standard pro- 
cedures with the MIDAS software: bias subtraction, flat- 
fielding with flats made on the dome or the sky, and the re- 
moval of cosmic-rays hits. The individual frames were then 



1 i \+ \ i j 


Tul 
j_ei. 


\-i Q Tl r\ A T H VTiAOIll'lO 

JDctllLl o6 H/XJJUbLlIL! 


Aff a 






(min) 


t"\ 
\ ) 


01037+1219 


VLT 


U (4) B (2) V (2) 


1.1 (V) 




HST 


F8UW (11.3) 


0.12 (F606W) 


09452+1330 


VLT 


V (120) 


0.7 (V) 


18240+2326 


OHP 


V (300) 


2.5 (V) 


18467-4802 


VLT 


U (8), B (2), V (2) 


1.1 (V) 




HST 


F606W (11.6) 


0.12 (F606W) 


20077-0625 


DT 


V (120) 


2.3 (V) 




HST 


F606W (11.6) 


0.12 (F606W) 


23320+4316 


OHP 


B (150), V (60) 


2.1 (B) 



a FWHM of field stars in the images 



shifted and stacked. In some cases, we corrected residual 
gradients in the background level using division by a nor- 
malized polynomial of order 1 or 2. 

The photometric calibration of the EFOSC2 observa- 
tions was made using photometric sequences close to the 
star T Phe and in the field SA 98-642 (Snodgrass HOTS)- 
For these observations, the relation between magnitudes 
and ADU counts is found to be: 

B = 25.80 - 0.20AT - 2.51og(F/t) 

for the B band, and: 

V = 25.86 - 0.11AT - 2.5 log(F/<) 

for the V-band, where X is the air-mass, and t is the ex- 
posure in seconds. For each relation, B or V is the stellar 
magnitude if F is the spatially integrated ADU count for 
the star, or the surface brightness if F is the ADU count per 
arcsec 2 measured on an extended object. The coefficients 
of X are the usual extinction coefficients for photometric 
nights at La Silla: the quoted values are from the documen- 
tation of the 1.54 Danish telescope. Overall the calibrations 
are expected to have an uncertainty of 0.1 mag, and the 
results are found to be consistent within 0.05 mag when 
comparing the two EFOSC2 runs. 

The photometric calibration for the EMMI observa- 
tions was made by considering un-saturated stars and their 
B-band magnitudes taken from the USNO-B1.0 catalog 
(Monet et al. 2005]); in this case the photometric uncer- 
tainty is ~ 0.4 mag (1-er). Comparisons of the calibration 
with the USNO catalog and the EFOSC2 photometric cal- 
ibration are consistent with this. 

The envelopes illuminated by the ambient Galactic radi- 
ation field are faint, typically a few percent of the sky back- 
ground. Hence the limiting factor for envelope detection is 
the sky noise. The exposures reported here are relatively 
long, and all the runs were made in dark time. The results 
are therefore representative of the quality of the data that 
can be obtained with this technique, using an intermediate 
sized, ground-based telescope. 

2.3. Additional observations 

For completeness, we include in the analysis the observa- 
tions from our earlier report (MH06) which are of com- 
parable quality to the NTT data. These observations are 
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Table 4. Observations of the circumstellar envelopes 



IRAS 


Tel. a 


Detection 6 


Envelope 
(B/V) 


S.B. C 

(Bo) 


Ellipticity 
(a/6) 


P.A. 

n 


a 

(") 


^cxp 

(yr) 


Core d 


00193-4033 


NTT 






<27.7 


<27.6 












01037+1219 


VLT 


St 


env 


24.3 e 


25.2 


1.07 + 0.03 


89 


5.2 


840 


pc 


04575+1251 


NTT 






<27.2 


<25.5 












05418-3224 


NTT 


St 




<26.9 / 


<26.8 












06012+0726 


NTT 






<25.6 e 


<24.7 












06255-4928 


NTT 


St 




<26.8 / 


<26.6 












07454-7112 


NTT 


st 


env 


25.8 


25.1 


1.12 + 0.03 


5 


3.9 


1290 


pi- 


09116-2439 


NTT 




env 


26.7 


25.9 


1.11 + 0.04 


47 


4.6 


1600 




09452+1330 


VLT 


St 


env 


25.0 9 


25.6 


1.03 + 0.02 


150 


45.8 


1770 


pe 


09521-7508 


NTT 


st 




<26.6 / 


<25.3 












10323-4611 


NTT 


st 


env 


25.9 


25.0 


1.20 + 0.04 


12 


4.1 


550 


pc 


12384-4536 


NTT 


st 


env 


25.8 


25.3 


1.14 + 0.03 


38 


3.6 


2320 




14591-4438 


NTT 




env 


25.5 


24.8 


1.10 + 0.03 


124 


2.7 


420 




16029-3041 


NTT 




env 


24.6 9 


23.8 


1.07 + 0.02 


41 


8.7 


4660 




17319-6234 


NTT 




env 


24.2 9 


23.9 


1.05 + 0.02 


65 


8.5 


2760 




18240+2326 


OHP 




env 


25.3 e 


25.6 


1.14 + 0.03 


46 


5.0 


1440 




18467-4802 


VLT' 1 




env 


23.7 e 


24.5 


1.03 + 0.02 


113 


5.0 


2310 




19178-2620 


NTT 


st 


env 


24.5 


24.0 


1.29 + 0.03 


156 


5.4 


3130 


pc 


20042-4241 


NTT 


st 




<27.6 / 


<27.4 












20077-0625 


DT 




env 


23.6 e 


24.3 


1.28 + 0.02 


58 


14.7 


3310 


pc 


23257+1038 


NTT 




env 


25.4 e 


25.9 


1.01 + 0.03 


157 


5.2 


3780 


pc 


23320+4316 


OHP 


st 


env 


26.8 


26.2 


1.13 + 0.06 


140 


7.0 


1760 


pc 



a Telescope used for photometry and imaging shown in Figs. 1-6 and Fig. A.l 
6 Detections in V or B: st = star, env = envelope 

c Surface brightness in mag. arsec -2 in B or V (designated in footnote) and effective Bo (see text) 
d Core: pc = polar core 
e V-band measurement 

^ Nominal limit based on extended source because of bright central star 
3 V surface brightness also measured: see section 3.3 
h VLT used for photometry, HST/ACS used for imaging 



summarized in Table 3, and have been processed in the 
same way. 

The observations include short exposures with the Very 
Large Telescope (VLT), longer exposures with the 1.54 m 
ESO Danish Telescope (DT) and the 1.20 m telescope of 
the Observatoire de Haute-Provence (OHP), and archival 
images made with the HST. The VLT-FORS raw images 
and calibration files were extracted from the ESO archive. 
The l^-band image of IRC+10216 was first shown by Leao 
et al. p006p . The HST data are from programs 9463 and 
10185, PI: R. Sahai. Further details of the observations are 
given in MH06. Observations in MH06 but not in Table 3 
include several non-detections which are much poorer than 
the NTT data, and AFGL 3068. In AFGL 3068, the illumi- 
nation is highly directional and the envelope shows a strong 
spiral pattern, which both limit the ability to measure the 
overall shape. 

3. The images 

3.1. Presentation of the images 

The main results of the observations are summarized in 
Table 4. Column (2) gives the telescope used for the V and 
B observations, and column (3) indicates the detection of 
the envelope and/or the central star (or the star-illuminated 
core) at these wavebands. For the NTT observations, seven 



of the envelopes are not detected and are discussed in the 
appendix. 

Figs. 1-4 present observations of each detected enve- 
lope, in three panels: a gray-scale image of the field in the 
waveband that best shows the envelope shape, a close-up 
contour map of the envelope, and the azimuthally averaged, 
radial intensity profile. 

In many of the gray-scale images it can be seen that 
field stars and galaxies lie along the lines of sight to the 
envelopes. In order to make the contour maps, the images 
are cleaned by removing these features using the follow- 
ing procedures. For bright stars whose halos extend over 
part of the envelope, we fit and subtract circularly sym- 
metric halo models. The models are derived from isolated 
field stars in the same image, and the fit is made by mini- 
mizing residuals in regions outside the envelope. The cores 
of these bright stars are also masked. In the case of IRAS 
07454—7112, a relatively bright elliptical galaxy lies close 
to the AGB star to the South-East, and this was success- 
fully subtracted using a Sersic galaxy profile (Sersic 1963). 
Dimmer stars and small galaxies, are simply masked. All 
masked areas are then replaced with local averages. This lo- 
cal average is determined by considering sectors centered on 
the AGB star with radial widths of 1-2", and typical ranges 
in position angle of ±20-30°. In this way, the masked areas 
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Fig. 1. Images and radial profiles of the AGB envelopes. Each row shows from left to right: the image in grayscale; a 
close-up contour map of the envelope, with field stars removed and their cores masked; and the radial profile of the 
envelope with the PSF. The field size is indicated in the images. The filters and contours, specified by the lowest contour 
(and contour interval) in units of the peak surface brightness of the envelope, are as follows (from top to bottom): IRAS 
01037+1219, V-band, contours 0.33 (0.43); IRAS 07454-7112, B-band, contours 0.31 (0.23); IRAS 09116-2439, B-band, 
contours 0.28 (0.26); IRAS 09452+1330, V-band, contours 0.32 (0.16). The images are square, with north to the top and 
east to the left. 




Fig. 2. Continuation of Fig. 1. From top to bottom: IRAS 10323-4611, S-band, contours 0.13 (0.25); IRAS 12384-4536, 
V-band, contours 0.25 (0.25); IRAS 14591-4438, 5-band, contours 0.35 (0.33); IRAS 16029-3041, V"-band, contours 
0.17 (0.17). The images are square, with north to the top and east to the left. 




Fig. 3. Continuation of Fig. 1. From top to bottom: IRAS 17319-6234, y-band, contours 0.12 (0.19); IRAS 18240+2326, 
V-b&nd, contours 0.34 (0.23); IRAS 18467-4802, F606W (HST/ACS), contours 0.27 (0.25); IRAS 19178-2620, S-band, 
contours 0.18 (0.22). The images are square, with north to the top and east to the left. 




Fig. 4. Continuation of Fig. 1. From top to bottom: IRAS 20077-0625, V-band, contours 0.20 (0.18); IRAS 23257+1038, 
V-band, contours 0.20 (0.31); IRAS 23320+4316, B-band, contours 0.29 (0.29). The images are square, with north to 
the top and east to the left. 



are filled with representative local intensities, and Gaussian 
noise similar to the sky noise is added in these areas. 

For the contour maps in Figs. 1-4, the cleaned images 
were lightly smoothed with a Gaussian of width 2 pixels 
(FWHM) with two exceptions noted in Sect. 3.5. The gray 
areas in the figures show the masked regions as described 
above. The solid contours are linearly spaced, and are speci- 
fied in the figure captions in units of the peak surface bright- 
ness of the envelopes; the equivalent values in mag arcsec -2 
are given in Table 4. In cases with a strong stellar core, ad- 
ditional dashed contours are included with a step size equal 
to the level of the highest solid contour. In Figs. 5 and 6 
(on-line only), the cleaned images are shown in color, with 
the contours superposed. 



3.2. General properties 

As described in MH06, deep optical imaging of AGB stars 
may reveal one or more of the following: the central star 
or an unresolved stellar core; the star-illuminated inner en- 
velope; and the extended envelope illuminated by the am- 
bient Galactic radiation field. The typical situation for the 
thick circumstellar envelopes is that the star and the star- 
illuminated core are faint (or not seen) in the B or V bands 
(because of extinction by the envelope), but the star be- 
comes dominant at longer wavelengths because of the de- 
creasing opacity of the dust envelope and the red spectrum 
of the star. 

For the star-illuminated inner envelope, the intensity of 
scattered starlight decreases rapidly with angular distance 



N. Mauron et al.: Imaging AGB circumstellar envelopes 



9 



from the center (~ 8~ 3 , Martin & Rogers Il987p . For exter- 
nal illumination, a faint, extended nebula is seen, with a 
shallow radial dependence. The externally illuminated neb- 
ula may also form a plateau of relatively constant intensity 
in the central regions where the optical depth to external 
radiation is £ 1. Examples of this (e.g., IRAS 16029-3041) 
may be seen in the radial profiles in Figs. 1-4. 

3.3. Surface brightness 

Column (4) of Table 4 lists the peak surface brightness in 
B or V for each circumstellar envelope, based on the pho- 
tometry described in Sect. 2.2 and MH06. For cases where 
the envelope is detected and the star is of low intensity 
or not detected, or the envelope plateau is extended, the 
measurement is straightforward. For cases where the enve- 
lope is detected, but the central star is of high intensity, 
we estimate the peak intensity by fitting and removing the 
stellar contribution at the center by using the profiles of 
field stars. For cases where neither the envelope nor the 
star is detected, we give upper limits assuming the source 
size is comparable to the seeing disk (FWHM). For cases 
where the star (or stellar core) is detected but not the enve- 
lope, a small source is essentially invisible; for these we give 
nominal limits based on the smallest source that would be 
detectable: 3" for IRAS 06255-4928 and 09521-7508 and 
5" for IRAS 05418-3224 and 20042-4241. 

For three envelopes we are able to measure both the 
B and V surface brightness. In addition to the B values 
given in Table 4, V = 25.04 mag arcsec -2 for the C-rich 
envelope IRAS 09452+1330, and V = 23.43 and 23.07 mag 
arcsec -2 , for the O-rich envelopes IRAS 16029-3041 and 
17319-6234, respectively. 

In addition to the measured values, Table 4 also lists an 
effective surface brightness in the S-band, Bo, corrected for 
interstellar extinction, and where necessary, transformed 
from U-band measurements. The extinction used is from 
the NErJ3 calculator, assuming all extinction along each 
line of sight lies between us and the star. For envelopes 
with only F-band measurements, we first correct for ex- 
tinction, and then adopt the extinction corrected color 
(B — V)o appropriate to the chemistry of the envelope. 
The colors are obtained from the three envelopes mea- 
sured in both bands given above. For the C-rich envelope, 
(B — V)o = 0.72 mag, and for the mean of the O-rich en- 
velopes, {B — V)o = 0.99 mag. Thus the two types of en- 
velopes are fairly similar. 

3.4. Envelope shapes 

Inspection of the images in Figs. 1-6 shows that the en- 
velopes exhibit relatively smooth and symmetric shapes 
that reveal their large scale geometry. They range from 
close to circular to distinctly elliptical. 

There are two additional features. First, three envelopes 
show elongated central components. These we interpret 
as the result of the leakage of light from the central star 
through polar (or bipolar) openings in the core of the en- 
velope, as is the case in IRC+10216 (Mauron & Huggins 
2000). Several additional cases are seen on smaller scales in 
HST images of the sample (MH06; Sahai l2007j) . The polar 
cores are noted in Table 4, and details are given in Sect. 3.5. 

http:/ /ned. ipac.caltech.edu 



A second feature of the envelopes of IRAS 16029—3041, 
17319-6234 and 18467-4802 is a one-sided asymmetry in 
the optically thick central region, which results from asym- 
metry in the external illumination. This is most clearly seen 
as a crescent shape in the highest level contours. In all three 
cases the direction of maximum illumination indicated by 
the crescent is within 30° of the direction of the Galactic 
plane. The asymmetry is reduced in the outer, optically 
thin regions of the envelopes where the observed intensity 
is less sensitive to the directionality of the illumination. 

In order to quantify the overall shapes of the envelopes, 
we fit ellipses to intensity contours at specified levels. Our 
procedure samples an intensity contour at intervals of 10° 
in position angle, omitting regions masked by the presence 
of field stars, as described in Sect. 3.1. An ellipse is then 
fit to the data using a non-linear least squares model, with 
the lengths of the semi-major and semi-minor axes (a and 
b), the center of the ellipse, and the position angle (P. A.) of 
the major axis as parameters of the fit. The shape is then 
given by the ellipticity E, where E = a/b. 

An example of the results of this procedure for the case 
of IRAS 17319-6234 is shown in Fig. 7. The figure shows 
data points sampled at levels of 20% and 40% of the enve- 
lope peak intensity, and the fitted ellipses for which values 
of E are 1.05±0.02 and 1.09±0.03, respectively. In spite of 
gaps in the data because of the masking of field stars in the 
image, the shapes are well constrained by the procedure. In 
this case the envelope is marginally non-circular, probably 
on account of asymmetric illumination seen in the plateau 
in the image in Fig. 5. 

Using this technique, we characterize the shapes of the 
envelopes on the largest size scales allowed by the observa- 
tions in order to minimize the effects of contamination by 
the central star (when visible) and the core features noted 
above. The results are given in Table 4. For cases with good 
S/N and resolved plateaus, we give the shape measured at 
20% of the peak. For the other envelopes, which are fainter 
and less well defined at this level, we give the shape at the 
40% of the peak; Table 4 lists the fitted values of E (±ler), 
a, and the P.A. of the major axis. The expansion time scale, 
^cxp corresponding to a is also given, based on the distances 
and expansion velocities given in Table 1. The envelopes are 
detected to significantly larger radii in the azimuthally av- 
eraged profiles, as seen in Figs. 1-4. 

The envelope shapes and orientations at higher intensity 
levels are generally similar to those in Table 4. However 
there are several exceptions where the inner ellipticity is 
larger and at a different PA, especially in cases which have 
the core features noted above. Individual cases are discussed 
below. 

3.5. Notes on individual objects 

IRAS 01037+1219 (IRC+10011). The extended envelope is 
essentially circular in shape. At the center a small scale 
polar core is seen in the HST/ACS F816W image shown 
in Fig. 1 of MH06, ~ 0'.'4 in length at P.A. -45°. There 
is also an asymmetry seen in CO in the central regions at 
approximately the same P.A. (Castro-Carrizo et al. 2010). 

IRAS 07454-7112 (AFGL 4078). The envelope at low 
intensities is approximately circular, but the shape is un- 
certain because of interference by field stars. Nearer the 
center, the shape is more definitely elliptical at P.A. +30°. 
The HST/ACS ^6061^ image of the star-illuminated cen- 
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ter (Sahai 2007) shows a polar core ~ 0'.'6 in extent at P.A. 

— 160°, which aligns with the central ellipticity. 

IRAS 09116-2439 (AFGL 5254). The envelope is con- 
sistent with circular, but the detailed shape is uncertain 
because of the large masked area of the field star near the 
center. 

IRAS 09452+1330 (IRC+10216). In spite of the high 
S/N of the image, the low intensity levels are affected by 
a gradient caused by a very bright field star outside the 
image shown (see Mauron & Huggins .1999). The shape in 
Table 4 refers to 30% of the peak intensity. At higher levels 
it is affected by shell structures which are slightly elliptical, 
and in the center there is a polar core at P.A. +12°, illumi- 
nated by the star (Skinner et al. 119981 Mauron & Huggins 
2000 ). Because of the large field, the contours in Fig. 1 are 
smoothed to a resolution of l'/6. 

IRAS 10323-4611. The envelope is elliptical at low 
intensities but more circular towards the center. The 
HST/ACS F606W image of the star (Sahai 120071) shows 
a polar core of ~ 0'.'4 at P.A. +25°. 

IRAS 12384-4536. Moderate but significant ellipticity 
at all levels. 

IRAS 14591-4438. Low S/N data. Roughly circular at 
low intensities, but may be elliptical in the center. 

IRAS 16029-3041 (AFGL 1822). High S/N data. The 
plateau shows asymmetric illumination, and a central dip 
characteristic of an optically thick core. The small but sig- 
nificant ellipticity in the extended envelope is likely caused 
by the asymmetric illumination. 

IRAS 17319-6234 (OH 329.8-15.8). High S/N data. 
The extended envelope is circular. The plateau shows asym- 
metry in illumination, and a central dip characteristic of an 
optically thick core. 

IRAS 18240+2326 (AFGL 2155). Small but uncertain 
ellipticity on account of bright, nearby field stars. 

IRAS 18467-4802 (OH 348.2-19.7). HST/ACS F606W 
image. The image and contours in Fig. 4 are smoothed to a 
resolution of 0'.'25. The extended envelope is circular. The 
plateau shows asymmetry in illumination, and a central dip 
characteristic of an optically thick core. 

IRAS 19178-2620 (AFGL 2370). Highly elliptical en- 
velope. Strong polar core at P.A. —155°, inclined to major 
axis at an angle of ^50°. 

IRAS 20077-0625 (AFGL 2514). Highly elliptical en- 
velope at all intensity levels. Fig. 8 shows a polar core at 
the center in the HST/ACS F606W image at P.A. -40°, 
roughly orthogonal to the major axis. 

IRAS 23257+1038 (AFGL 3099). The extended enve- 
lope is circular and centered on the AGB star. There is a 
strong polar feature at the center, at P.A. —175°. 

IRAS 23320+4316 (AFGL 3116). Low S/N data. The 
envelope is roughly circular, with evidence for an elliptical 
core at P. A. +1 10°. The HST/ACS F606W image of the 
star (Sahai [2007P shows a polar core ~ 0'/3 in extent at P.A. 

— 150°, approximately orthogonal to the elliptical core. 

4. Discussion 

4.1. Envelope surface brightness 
4.1.1. Detectability of the envelopes 

The detectability of an envelope in scattered ambient light 
depends on the effective optical depth for scattering, the 




I i i i i i I 

400 440 480 

AX (pixels) 

Fig. 7. Contour fitting for IRAS 17319-6234. The inner 
(grey) and outer (black) crosses sample the envelope at 10° 
intervals in azimuth and at 20% and 40% of the peak inten- 
sity, respectively. Gaps in the data correspond to masked 
regions of the image. The curves show the fitted ellipses; 
the small inner circles show the fitted centers. 




Fig. 8. Polar core at the center of IRAS 20077-0625. HST 
F606W image. The data have been binned to a resolution 
of O'.'l. Field size 4"3 x 4'.'3. 



strength of the radiation field, and the interstellar extinc- 
tion along the line of sight. MH06 describe how the envelope 
opacity is conveniently parameterized by an estimate of the 
angular size of the optically thick plateau region using: 

AG = 5" ( M )( - d )~ 1 

V 10- 5 M Q yr- in 10km s- i; V 1000pc y ' 

If this exceeds the resolution of the observations, the center 
of the envelope should attain maximum surface brightness. 
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Fig. 9. Surface brightness vs. galactocentric radius. Filled 
and open symbols denote C-rich and O-rich envelopes, re- 
spectively. Circles are detections, triangles upper limits on 
the surface brightness. 

Based on the parameters listed in Table 1, all the en- 
velopes observed have AO 0'.'8, as given by the above 
equation. They should therefore all be partly or fully re- 
solved. If we consider the envelopes with good measure- 
ments (ignoring 4 with nominal upper limits due to the 
presence of a bright central star), there are 15 detections 
and three upper limits. One of the limits is sensitive, but it 
belongs to the envelope with the lowest A9, so it it proba- 
bly just too small to be seen. The two other non-detections 
are for the two objects with the highest interstellar extinc- 
tion (1.6 and 2.0 mag. in B) which prevents reaching low 
intensity levels. Hence, if allowance is made for extinction, 
the characterization in terms of A9 is a successful guide to 
detectability. 

4.1.2. Envelope surface brightness and the interstellar 
radiation field 

After correction for extinction, the peak surface brightness 
Bo shows a range of more than 2 mag (see Table 4). The 
scatter does not correlate with latitude (b) or height above 
the Galactic plane (z), and is not clearly related to A9 or 
M. However, it does appear to depend partly on envelope 
chemistry and longitude (I), and both are found to result 
from a remarkable correlation with Galactic radius R. 

Fig. 9 shows the correlation of Bq with Galactic radius, 
with R calculated using the stellar distances and coordi- 
nates given in Table 1, and assuming R = 8.0 kpc for 
the Sun. It can be seen that the surface brightness de- 
creases rapidly with R with an approximately linear re- 
lation. Omitting the upper limits for reasons mentioned 
above, a least squares fit gives: 

B = 0.90 (±0.17) R+ 18.20 (±1.26) 

and a residual rms scatter in Bq of 0.41, in units of 
mag arcsec -2 . We interpret this relation to be the effect of 
a radial gradient in the interstellar radiation field. The lin- 
ear dependence in magnitudes translates to an exponential 



law in the mean interstellar intensity. As far as we know, 
this is the most direct measurement so far of the gradient 
in the solar neighborhood. Our data suggest that the ra- 
diation field decreases by a factor of ~5 between R = 6 
and R — 8 kpc. This is fairly consistent with the model of 
Mathis et al. (1983) who find a decrease of a factor of 7 
between R = 5 kpc and the solar radius. A smaller factor 
of 3 is favored by the more recent study of Sodroski et al. 
(1997). 

Inspection of Fig. 9 shows that the C-rich envelopes 
(filled symbols) follow essentially the same relation as the 
O-rich envelopes (open symbols). It would appear that the 
scattering properties of these different envelopes in the B- 
band are rather similar, and this is consistent with the sim- 
ilarities in the B — V colors measured in Sect. 3.3. We 
therefore interpret the observed difference in the average 
brightness of O-rich and C-rich envelopes in terms of their 
locations in the Galaxy. In our sample of AGB stars, the 
C-rich envelopes lie on average at larger values of R, and 
this is consistent with the higher relative fraction of carbon 
stars in the outer Galaxy (e.g., Jura et al. I1989|) . At larger 
Galactic radii, the C-rich envelopes are illuminated by a 
less intense radiation field. 

4.2. Envelope morphology 

4.2.1. Distribution of observed shapes 

All the envelopes in Figs. 1-4 show relatively smooth radial 
profiles except for IRC+10216, which has sub-structure in 
the form of multiple arcs. However, IRC±10216 is much 
nearer than the other objects (see Table 1), so these could 
have similar structure which is not resolved at the larger 
distances. None of the envelopes show evidence for major, 
single-shell ejections which are seen in a few AGB envelopes 
and are thought to be caused by helium shell flashes (e.g., 
Olofsson et al. I^U0D|) . 

As described in Sect. 3, the large scale shapes of the 
envelopes are roughly circular or elliptical, consistent with 
spherically or axially symmetric three dimensional struc- 
tures seen in projection. The observed shapes are con- 
veniently characterized by the ellipticities (E) listed in 
Table 4; recall, E = a/b (> 1.0), where a and b are the semi- 
major and semi-minor axes, respectively. The distribution 
of ellipticities for the sample (in order of increasing E) is 
shown in Fig. 10, which also indicates the formal l-a and 
3-<7 uncertainties. As explained in Sect. 3.4, the ellipticities 
we use to characterize the distribution refer to the over- 
all shapes of the extended envelopes. If we use the shapes 
measured closer to the star (e.g., at 60% of the peak inten- 
sity), the distribution is similar to that in Fig. 10, but with 
somewhat higher ellipticities at the top end due mainly to 
the contribution of the core features in a few cases. 

From Fig. 10, it can be seen that deviations of the over- 
all shapes of the envelopes from circular symmetry are rel- 
atively small. Approximately half of the sample (47%) have 
E 1.1; 20% have E 1.2, and the rest are in between. 
The observed shapes of axially symmetric envelopes depend 
of course on the inclination angles, and they will appear 
circular when viewed end-on. However, there is no known 
bias in the inclinations of the observed envelopes, so we can 
assume them to be randomly oriented, with a median in- 
clination of 60° to the line of sight, and rarely seen end-on. 
On this basis we can draw the following general conclusions 
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Fig. 10. Distribution of the measured large scale elliptici- 
ties of the AGB envelopes. The vertical error bars denote 
the formal 1-a and 3-cr uncertainties. 

on the 3-dimensional shapes of AGB envelopes: 1. A signifi- 
cant fraction (~ half) are close to spherically symmetric. 2. 
At least 20% are significantly flattened. 3. Highly flattened, 
disk- like envelopes are relatively rare or non-existent. The 
results of Neri et al. (| 1998[) based on observations of AGB 
envelopes in the millimeter lines of CO are in qualitative 
accord with our conclusions. 

4.2.2. Shaping mechanisms 

There are two separate categories of shaping that can af- 
fect the morphology of the circumstellar envelopes: intrinsic 
shaping by or close to the central star, and interaction with 
the ISM. The latter can produce significant effects on large 
angular scales, typified by a one-sided bow shock due to the 
relative motion of the star and the surrounding medium. 

We first consider whether the ISM interaction is likely 
to be important for the observed shapes by comparing es- 
timates of the bow shock stand-off distance Rq for each 
object, taking into account the height above the galactic 
plane, using equations (2) and (5) of Cox et al. ((2,012), the 
stellar data in Table 1, and assuming a typical velocity rel- 
ative to the ISM of 30 km s _1 (Feast & Whitelock I20"0"0"f . 
We find that i?o for the envelopes observed ranges from ~ 
10 18 to 10 20 cm, and exceeds the size scale corresponding 
to the angular radii (a) given in Table 4 by a median factor 
of 70, with a range of 30 to 3000. We conclude that interac- 
tion with ISM is unlikely to have a significant influence on 
the measured shapes. The observed morphologies are also 
consistent with this conclusion. 

Our understanding of the intrinsic shaping mechanism 
of AGB envelopes is not well developed, but the expectation 
is that envelope shapes are different for single stars and 
binaries. Companions to AGB stars are, however, extremely 
difficult to detect (Jorissen & Frankowski I2008|) and it is 
probable that hidden binaries are included in our sample. 

For single AGB stars, which rotate very slowly, the 
zeroth order expectation is for spherically symmetric en- 
velopes, although the extent of deviations caused by mag- 



netic fields or related phenomena is not known. For AGB 
stars in binaries, the companion is expected to influence 
the shape in several ways. Close-in companions can cause 
tidal effects and spin-up, and modify the effective gravity of 
the star (e.g., Frankowski & Tvlenda l200ip . However, two 
additional effects are likely to be more important overall 
because they operate over a much wider range of separa- 
tions: the flattening of the envelope by the gravitational 
field of the companion, and the effects of jets powered by 
accretion onto the companion. 

Gravitational flattening. The gravitational flattening of an 
envelope towards the orbital plane of a companion is a 
prominent feature of gas dynamical simulations of binary 
AGB systems (Mastrodemos & Morris [1999} Gawryszczak 
et al. I2002[) . Low mass companions or wide separations lead 
to nearly spherical envelopes, and close-in, high mass com- 
panions lead to flattened envelopes. 

As a rough guide to the relation between envelope shape 
and binary parameters, Huggins et al. (|2009|) have obtained 
an approximate formula for the flattening in terms the com- 
panion mass M c , the wind speed V, and the binary separa- 
tion d, based on the simulations of Mastrodemos & Morris 
(|1999p . They use a highly simplified, axially symmetric 
model envelope with a density ~ r~ 2 , and characterized 
by a shape parameter K, where K — n e (r)/n p (r) is the 
ratio of equatorial to polar density; the relation between K 
and the binary parameters is given in their equ. (4) . For an 
envelope illuminated by external radiation, they also give 
a prescription for the observed shape for a specified incli- 
nation of the symmetry axis to the line of sight (i). The 
observed shapes are given in terms of the ratio of column 
densities on the major and minor axes, N ma j(9)/N min (0), 
but this can be shown to be equivalent the ellipticity (E) 
used in this paper. 

To illustrate the range of binaries that could give rise 
to the observed shapes, Fig. 11 shows contour maps of syn- 
thetic images based on the simplified envelope model. The 
left hand panel shows the dependence of the observed shape 
on M c and i for d = 10 AU. The right hand panel shows 
the dependence on M c and d, for the median inclination 
(i = 60°). A representative value of V = 10 km s^ 1 is used 
for all cases. 

It can be seen from the figures that the model produces 
envelope shapes that are qualitatively similar to those ob- 
served, i.e., with round or approximately elliptical contours. 
(The pinching effect on the minor axis of the flattest model 
is an artifact of the specific latitude dependence of density 
in the simple model, and is not expected to be a feature 
of real envelopes.) For masses ^ 0.1 M Q , the companion 
has essentially no effect at any angle, and the envelopes are 
very close to circular. For typical inclinations (i ~ 60°), 
companions with masses from 0.2 to 1.0 M©, significantly 
affect the envelopes, producing E 1.2 for separations up 
to 10 and 40 AU, respectively. Thus the range of observed 
envelope shapes is consistent with a spectrum of AGB stars 
ranging from single stars to binaries with relatively close 
companions of moderate mass. 

Polar cores and jets. The idea that binary companions 
are responsible for shaping the non-spherical envelopes is 
strengthened by the presence of the polar core asymme- 
tries noted in Sect. 3. The fact that these elongated regions 
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of enhanced radiation are seen in the central regions of ap- 
proximately half of the sample (Table 4) is striking. 

The best studied case is IRC+10216 (IRAS 
09452+1330) where the effect is caused by the leak- 
age of stellar radiation along bipolar axes close to the star, 
with one side much brighter than the other because of 
extinction (Skinner et al. UME[ Mauro n fe H uggins |2"0"0"0"|) . 
Other cases are probably similar fSahai l2"0TJTj) . although in 
the absence of spectroscopy, localized emission cannot be 
ruled out. 

The polar cores are not expected to be features of the 
mass loss of single stars, but they are expected in binary 
systems. One production mechanism is the diminution of 
mass loss along the polar directions which is an integral 
part of the binary shaping process, as seen in the simula- 
tions of Mastrodemos & Morris (|1999p . A second, related 
mechanism is the sculpting effect of jets from an accretion 
disk around the companion (e.g., Soker & Rappaport 2000). 
The effects of jets are seen in the Mira system (Josselin et 
al. 120001) . and in a more extreme form in proto-PNe (e.g., 
Huggins et al. 12004)1 . These mechanisms can act in parallel, 
and the dominant effect could vary from object to object. 

Since polar cores and gravitational shaping likely arise 
from similar circumstances (i.e., the presence of a com- 
panion) there is good reason to suppose that they could 
be correlated. Our observations provide evidence for this 
view. The three envelopes with the highest cllipticitics all 
have polar cores (Table 4). The polar cores are present in 
about half the sample, somewhat more than the fraction 
that are clearly flattened. It could be that the polar cores 
are formed by low as well as moderate mass companions, 
and are therefore more widely detectable. The relation be- 
tween these features clearly warrants further study. 

The high frequency of occurrence of polar cores in 
the observed sample also places interesting constraints on 
their physical origin. If they are intrinsic to envelope shap- 
ing, they will be permanent features of the envelopes. 
However, if they are caused by jets, the time evolution needs 



to be considered. Several cases where jets are suspected 
are discussed by Sahai (|2"P0T)) . and IRAS 010 37+12 19 
(IRC+10011) has been modeled by Vinkovic et al. (j2HM|) in 
terms of jets which have very recently turned on (< 200 yr) 
and are in the process of sculpting the envelope. Most of 
the cores observed here are small in radial extent, and with 
the high speed of jets, have comparably short time scales. 
Since we observe polar cores to be common, it is statisti- 
cally unlikely that they have all formed within the last few 
hundred years. For the jet interpretation, the jets are lik- 
ley to be weak or intermittent, so that the effects are local 
to the star. In this picture the more powerful jets seen in 
proto-PNe would develop later, when the companion inter- 
acts more strongly with the AGB star and accretes at a 
higher rate (Huggins I2007|) . 

4.2.3. Population synthesis 

The question of binary shaping can be investigated further 
by comparing our observations of envelope shapes with the 
results of a recent population synthesis of AGB stars by 
Politano & Taam (|201 1|) . For a population consisting of 
single stars and binaries, they have calculated the distribu- 
tion of the gravitational flattening parameter K of Huggins 
et al. (|2009p described above. We have used their results 
to predict the distribution of the observed shapes of en- 
velopes, assuming random orientations. Any contribution 
to the shaping by jets would be additional. We report the 
results in Fig. 12, using the median values of E in 15 quan- 
tiles, to represent the sample of observations. 

The population synthesis assumes 50% of stars form as 
single stars and 50% from in binaries on the main sequence. 
For the binaries, we consider two representative cases for 
the assumed distribution of qo, the ratio of the secondary to 
primary mass on the main sequence. Case 1 is the standard 
case of Politano & Taam (|201ip with a flat mass distri- 
bution, g (go) = 1- Case 2 is a distribution that favors low 
mass companions with g(qo) — 9o~°' 9 - F° r case 1j 32% of the 
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Fig. 12. Comparison of the observed ellipticities of the 
AGB sample with calculated distribution functions. The 
points are the observations, and the vertical error bars 
denote the 1-a uncertainties. The upper and lower gray 
histograms are population synthesis predictions including 
the measurement errors for case 1 (g{qo) = 1) and case 2 
(5(90) = %°' 9 )> respectively. The upper and lower black 
histograms are population synthesis predictions including 
a residual ellipticity for case 1 and case 2, respectively. See 
text for details. 



envelopes on the AGB are perturbed by a companion, and 
68% are unperturbed (either single stars, or with distant 
companions); for case 2 these numbers are 18% and 82%, 
respectively. 

If single stars produce exactly spherical envelopes, both 
cases 1 and 2 would have more than half the AGB en- 
velopes with nominal ellipticities of E — 1.0 (in the absence 
of measurement errors). For realistic comparisons with the 
observations we include two effects. The gray histograms 
in Fig. 12 show the predicted distributions, where we in- 
clude by simulation the effect of the average uncertainty 
in fitting the ellipses a — ±0.03. The black histograms in 
Fig. 12 show the two cases if we assume a residual ellip- 
ticity, normally distributed with mean 1.06 and dispersion 
±0.06 (1-cr) for all envelopes, in quadrature with any bi- 
nary perturbation. This residual effect is chosen for illus- 
tration, and could arise from a combination of systematic 
effects such as residual asymmetric illumination and plate 
gradients, and as a natural dispersion in the true shapes of 
unperturbed envelopes. 

In spite of the extreme simplicity of the underlying enve- 
lope model and the complexity of the population synthesis, 
it can be seen in Fig. 12 that the synthesized curves exhibit 
the qualitative characteristics of the observed distributions. 
The simulations suggest that the envelopes with ellipticities 
close to 1.0 are mainly unperturbed envelopes with a small 
dispersion in their shapes; and the less frequent, larger el- 
lipticities are the results of binary interactions. There are 
significant quantitative differences in the predicted curves 
for the different assumptions about the binary mass ratios. 
It is encouraging that the observational data fall between 



the two cases shown, but the observational sample is too 
small to discriminate between them at this stage. 

From an empirical point of view, if the polar cores indi- 
cate binaries as we suggest, the observation of cores in half 
the sample would indicate a somewhat greater percentage 
of AGB binaries than predicted by the population synthe- 
sis. A larger sample would be important to investigate this 
further. 



4.3. AGB stars as precursors of PNe 

The evidence for a hidden population of AGB binaries pro- 
vided by our observations is of interest for understanding 
post- AGB evolution. In particular, the importance of bina- 
ries in shaping planetary nebulae (PNe) is a subject that 
is widely debated (De Marco I2009p . and it is uncertain 
whether a binary companion is an essential ingredient for 
PN formation. It is observed that 12-21% of PNe have close 
binary companions (Miszalski et al. I2009|) and these have 
evidently passed through a common envelope phase. Other 
companions at larger distances may be present, but have so 
far proved elusive to observation (e.g., Hrivnak et al. 1201 lj) . 

The objects in our sample of AGB stars with high mass 
loss rates are expected to be the immediate precursors of 
PNe. We find that ^20% of the envelopes already show 
significant shaping effects that we attribute to the pres- 
ence of substantial companions. We also find evidence for 
additional, more subtle shaping effects in the form of po- 
lar cores in up 50% of envelopes, which also probably arise 
from companions (including low mass companions). These 
numbers are preliminary because the sample size is small, 
but they demonstrate that detailed observations of AGB 
envelopes have the potential to characterize aspects of the 
interactions that contribute to the shaping of PNe. 



5. Conclusions 

We have carried out a deep imaging survey of 22 AGB stars 
with high mass loss rates in order to investigate the geom- 
etry of their circumstellar envelopes. We report the detec- 
tion of 15 envelopes in dust-scattered Galactic radiation, 
and we characterize their properties in terms of the surface 
brightness, radial intensity profile, and observed shape. 

We find that the surface brightness of the envelopes 
shows a rapid decrease with Galactic radius, which we in- 
terpret as a steep gradient in the interstellar radiation field. 
As far as we know, this is the most direct observation of 
this gradient in the solar neighborhood. 

The envelopes show a range of geometries: approxi- 
mately half are close to spherically symmetric, and ^20% 
are distinctly elliptical. We interpret the shapes in terms 
of populations of single stars and binaries whose envelopes 
are flattened by a companion. The observed distribution 
of the ellipticities of the envelopes is qualitatively consis- 
tent with the results of population synthesis models. We 
also find that approximately half the sample of envelopes 
exhibit small-scale, polar cores which we interpret as the 
escape of starlight through polar holes produced by com- 
panions. 

Our observations of envelope flattening and polar holes 
point to a hidden population of binary companions within 
the circumstellar envelopes of AGB stars. These compan- 
ions affect the geometry of the envelopes on the AGB, and 
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are potentially important guides to the evolutionary chan- Vinkovic, D., Blocker, T., Hofmann, K.-H., Elitzur, M., & Weigelt, G. 
nels that lead to binaries observed in the post- AGB phase. 2004, MNRAS, 352, 852 

Whitelock, P., Menzies, J., Feast, M., ct al. 1994, MNRAS, 267, 711 
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Appendix A: Envelope non-detections 

Images of the fields for which the AGB star but not 
the envelope are detected are shown in Fig. A-l. IRAS 
06012+0726 (AFGL 865) is not detected in V or I, and is 
not shown. It is however well detected in J, H, and K in the 
2MASS survey, with J = 13.38, H = 9.28, and K = 6.06. 
At a Galactic latitude b = — 7, this object has the largest 
extinction (Ay ~ 1.6 mag.) based on the NED extinction 
calculator, which probably accounts for the non-detection 
at short wavelengths. 
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Fig. A.l. Images of objects with non-detection of the envelope. Ticks indicate the AGB star. Top row, left to right: 
IRAS 00193-4033 (AFGL 5017), EFOSC2 ^-band image. IRAS 04575+1251 (AFGL 5134), EFOSC2 F-band image. 
IRAS 05418-3224, EFOSC2 B-band image. Bottom row, left to right: IRAS 06255-4928, EFOSC2 5-band image. IRAS 
09521-7508, EMMI B-band image. IRAS 20042-4241, EMMI S-band image. 
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Fig. 5. Images of the circumstellar envelopes (on-line only). The images are shown in color with the field stars replaced, 
and with the contours of Figs. 1-4 superposed. Top row, left to right: IRAS 01037+1219, I^-band; IRAS 07454-7112, 
B-band; IRAS 09116-2439, B-band. Middle row, left to right: IRAS 09452+1330, V-b&nd; IRAS 10323-4611, B-band; 
IRAS 12384-4536, 7-band. Bottom row, left to right: IRAS 14591-4438, B-band; IRAS 16029-3041, V^-band; IRAS 
17319—6234, 1^-band. All images are square, with north to the top and east to the left. 
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Fig. 6. Continuation of Fig. 5 (on-line only). Top row, left to right: IRAS 18240+2326, V-band; IRAS 18467-4802, 
F606IF; IRAS 19178-2620, B-band. Bottom row, left to right: IRAS 20077-0625, F-band; IRAS 23257+1038, F-band; 
IRAS 23320+4316, i?-band. All images are square, with north to the top and east to the left. 



